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SUMMARY. 


The TBLLUS Project entered its central application phase. 

A simple procedure to evaluate actual evaporation vas derived 
by linearizing the surface energy balance equation, using 
Taylor's expansion. The original multi-dimensional hyper- 
surface could be reduced to a linear relationship betveen 
evaporation and surface temperature or to a surface rela- 
tionship involving evaporation, surface temperature and al- 
bedo. This procedure permits a rapid sensivity analysis of the 
surface energy balance equation as well as a speedy mapping 
of evaporation from remotely sensed surface temperatures and 
albedo. Comparison with experimental data yielded promising 
results . 

Testing the validity of evapotranspiration and soil moisture 
models (as TERGRA) in semi-arid conditions is one of the aims 
of the TELLUS Project. It corresponds with the Commission's 
initiatives in favour of the Associated Developing Cotintries. 
In particulax' it can be viewed as an initial contribution to 
the solution of menacing desertification problems. Wheat was 
the crop chosen for a continuous measurement campaign made in 
the South of Italy. Radiometric, micrometeorologic, agronomic- 
and-soil data were collected and centralized at JRC for pro- 
cessing and interpretation.. 

A number of studies on the relations between HCMM sensed sur- 
face temperatures and different land-use pattern as well as 
local topography were performed. 


HCMM night digital data recorded over Belgiiun and correspon- 
ding LANDSAT data were analyzed. Forested areas in the Ardennes” 
could be veil delineated by their high surface temperatures. 

Prom an HCMM scene of the Upper Rhine Valley between Basel 
and Frankfurt it was confirmed that HCMM's resolution is not 
sufficient to determine the boundaries of various geographic 
units sm-^ler than 600 m or to establish the exact limits of 
such tf ' aphic structures. 

On the contrary, HCMM thermal images were found to be much 
more comprehensive for studies on a regional scale and ther- 
mal patterns at this scale can be related to and explained 
by topographic structure. HCMM's surface temperatures appear 
to be well correlated with forest distribution. Forests on 
steep slopes or convex terrain show very high temperatures 
(as warm as city centers). Contrast between forested and un- 
forested terrain is the most essential factor of the distri- 
bution of warm and cold surfaces. Many factors seem to in- 
fluence the distinctive behaviour of forested surfaces (pre- 
sence, inversion or interruption of boundary layers, sensible 
heat flux between air and forest being higher than adjoining 
grassland and arable lands, etc.). Presence of cold air re- 
servoirs and cold air streams was also explained for different 
situations. 

Mesoscale heat budgets particularly affected by anthropogenic 
heat release were examined. This investigation was performed 
on a segment of an HCMM scene encompassing the Upper Rhine 
Valley between Basel and Frankfurt and the surrounding highlands. 
The influence of urban settlements of various sizes was inter- 
preted as well as the presence of industrial complexes which 
were identified and analyzed in detail, 

A further calibration of the HCM radiometer was attempted 
using water tem.peratures of Lake Geneva and the radiation 
balance of a pine forest. These attempts cover the period of 
May-November 1978 and a temperature range from 9-23® C. 
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From their results one may tentatively conclude that the 
HCMR underestimated true surface temperature by about 5®C* 
Additional calibration studies are needed to reach a more 
definite conclusion. 

The RADTRA model for atmospheric corrections was employed 
extensively in the calibration study. The model had been mo- 
dified by the insertion of a correct value of the mass absorp- 
tion coefficient considerable difficulties were encoun- 

tered with the model. 

These problems were communicated to the HCMM Project Scientist 
and the model has been reprogrammed since. 

Some of the above-mentioned activity, already accomplished 
by the authors, is described in the annexed documents forming 
part in extenso of this report (see list at page 39) a 


1, INTRODUCTION 


TSLLUS is a remote sensing project centered arovmd the satellite 
EXPLORER-A, bearing the sensors HCMM (Heat Capacity Mapping Mis- 
sion) • The special feature of this satellite measuring the tempe- 
rature cycling (n/d) of objects on the Earth surface makes it 
suitable for hydrological and/or geological prospection. 

The main objective of TELLUS is to demonstrate that such tempera- 
ture cycling, or diurnal temperature only, can be related to 
moisture content in bare and vegetation-covered soils and to an- 
thropogenic heat release. The study of the influence of natural 
modifications on the regional h'^at budget is another important 
subject of investigation. 

The objectives of the proposed programme, vhile partially fulfil- 
ling a number of specific needs in the policy of agriculture and 
the environment, towards European and Associated African Countries 
also correspond to the specific research objectives of the parti- 
cipating Institutes. 

The role of the JRC in the TELLUS project consists of two distinct 
parts; 

- A coordinating activity for the contribution from national in- 
stitutes (co-investigators) working on different test sites and 
topics (see the 1st Progress Report). This coordination includes 
the present liaison with NASA, Space Goddard Flight Center. 
(GSFC) . 

- experimental activity, essential in itself, on test site 
No. 1, and orgcinisation and participation in Joint Flight Cam- 
paigns. In addition the JRC is dispatching, processing and in- 
terpreting data received from NASA/GSFC an.d BSA/EARTHNBT, Cen- 
tre de M4t4orologie Spatial de Lannion. 

2. TELLUS FRAMEWORK 

The main lines of TELLUS' organisation and objectives are here 
recalled. This will help to better understand who did, why and 
where dwing the reporting period. For a more detailed informa- 
tion please look it up in the 1st Progress Report (Chapt.s 2,3 
and 4) and in ref f ^ J and £ 2 J\ 
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The Co-investigators 

The names and abbreviations of the University Departments and In- 
stitutes participating in the TSLLUS mission are listed in Table 1 • 
Coinvestigators framevork consists of six groups. The JRC Ispra 
and some EC’s Directorates make up the Pi’s leading group. 

2.2. The Test-sites 

Geographic as well as technical considerations led to arrange in 
five groups the fifteen test-sites as originally proposed by the 
Coinvestigators. No particular test-site was proposed by the JRC 
Ispra within TELLUS. His research staff works in collaboration 
v/ith TELLUS coinvestigators on their national test-sites. A test- 
site coordinator was designated for each group according to the 
Institutes involved. 

Arrangement in groups of TELLUS Test-sites is reported in Table 2. 
Pig. 1 shows the distribution of the test-sites proposed by Coin- 
vestigators within the overall test-site "Europe'* which was agreed 
by JRC and NASA (boundary in brackets). 

2.3. Mission’s objectives 

According to the interest of the EC Commission and the research 
proposals made by the Coinvestigators, there are three main the- 
matic lines of investigation; 

1 • Evaluation of evapotranspiration and moisture content of bare 
soils and of soils covered by vegetation. 

2. Study of the interaction between natural phenomena and meso- 
scale heat budget. 

3. Man-made changes and their impact on regional heat budgets. 


Table 2 ; TSLLUS test sites, arranged in groups 


Test-Site Group 

Test Site 

ITALY 

Puglia and Basilicata* 

Sardegna 

Emilia 

FRANCE 

B ouches du Rh3ne 
Bretagne 

GERMANY 

Rhine Valley i 
Rhine Valley 2 
Northern A3.ps 
Northern Germany 

UNITED 

KINGDOM 

England i 
Vales 

England 2 
England 3 

BENELUX 

Benelux 

DENMARK 

Greenland** 


Basilicata proposed by UR/DG 
No HCMi’1 receiving station 






















3. DATA RECEPTION AND DISTRIBUTION 

JRC receives the HCMM data from two origins: 

1 ) from NASA/GSFC 

2) from CMS Lannion (BARTHNET) 

For more details please refer to the 1st PR, Chapt.5 pg.9 
and the 2nd PR, Chapt.3, pg.3. 

3.1 NASA/gsfc data delivery 

3.1.1. Standing Order data. Situation at Aug. 31 st, 1980 

At this moment of TELLUS' inves-igation no more Standing Order 
(so) data can be waited for at the JRC and we may draw an overall 
conclusion of the normal SO distribution. 

Historically, let us briefly recall the flow rate of SO data 
from Greenbelt to Ispra, 

In Autumn 78 preliminary data arrived first, on form of tran- 
sparencies and tapes, whereas the SO data (DVIS, DIR, NIR) 
started arriving early March 1979 on form of neg. transparencies 
(24 cm X 24 cm). 

The well known delay of Standing Order data led to the esta- 
blishment of the Priority Processing List (PPL) agreed by all 
Principal Investigator and NASA in Jan. 1979. 

In spite of such initiative the NASA response to the PPL was 
surprisingly poor. At the end of March 1980 nearly 56 only 
of the PPL was satisfied (see Table 6). 

Now, at the end of the final distribution the rejjponse to the 
PPL is increased slighty in respect to the preceding one (see 
Table 7). As a matter of fact it does still remain less than 
63 % although one must recognize that every effort has been 
made by the HCMM staff to compensate the delay. 

Up to now the JRC has received sets of scenes of 368 orbits, 
but the 38 processings of 370 scenes made by NASA from March 
till August did not satisfy exhaustivelly the PPL. 

It is therefore understandable that various TELLUS Co-Inve- 
stigators (such as the French ones) waited till the end of 
SO distribution -as formally established - to decide whether 
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TABLE 3. STATUS OF PRIORITY PROCESSING LIST AT MARCH 31 St. 1980 
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TABLE 4. STATUS OF PRIORITY PROCESSING LIST AT AUGUST 31 St, 1980 
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to keep to the PPL or try to look for other Night/Day pairs. 

Some others instead gave up the first one and repetitively 
asked for Retrospective Orders ( RO) • 

The Table 5 does reflect the actual distribution of Standing 
Order data made via Ispra. A comparison between the situation 
of the responses to the PPL and the actual RO distribution 
per country may offer an interesting indication. 

3.1.2 Retrospective Order Data 

Retrospective Order data were requested by many Coinvestiga- 
tors after reception of Standing Order data. 

CCT's of 101 scenes (DVIS-DIR-NIR) have been requested and 
95 of them have been received by the JRC and distributed. 

Four "scenes" with Temperature Difference (TD) data and 
Apparent Thermal Inertia (ATI) data have also been requested 
and received. 

The overall RO situation has been resumed in Table 6. 

3.1.3 Priority Coverage. 

A series of priority coverage has been requested by various 
Coinvestigators to match with concommittant ground measurements 
or to establish an interccmparison with other satellite data 
(see list of Table 7). 

The setting up has however revealed difficult owing to different 
misunderstanding. In 1978, it was difficult to establish whether 
the satellite could receive or not the data, because of conges- 
tion at the Madrid Station from other satellite transmissions. 
Sometimes no switching could be organised at all. 

After March 1980 lowering of orbit gave difficulty to establish 
what was the actual track, as explained in 3.1.4. 

3 9 1 .4 Lowering of satellite's orbit s 

Finally let us consider an important factor affecting the dis- 
tribution of the HCMM information: lowering of HGMM sat<='Xlite 
which disturbed the whole 16 tracks cycle. The exact modification 
and consequences of this was not clearly specified in time to Co- 
investigators when requesting Priority switch-on for special 
HCMM coverages to NASA. This inconvenient caused uncertainties and 
difficulties in organi:!.ation of related ground-truth activities. 
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Sibari 

24.07.78 

39.44/1 6, a?E 
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3.2 BSA/SARTHNBT data delivery 

Prom May 1979 CMS Lannion has been producing and distri- 

buting regularly Quick-Looks for an amount of 527 orbits. 

Table 3 shows the distribution of the Quick-Looks orbits by 
CMS Lannion to TELLUS Co-investigators. 

Several Co-investigators after reception of Quick-Looks, 
requested the Retrospective Order data using procedure of 
standard cutting of the orbits (see 2nd PR, Chapt. 3.2.2). 
Retrospective Orders were requested for both transparencies 
and CCTs, CCTs have been ordered for an amount of 74 scenes, 

71 of which have been received by JRC and distributed (see 
Table 6 ). 

The good organi7.ation of CMS Lannion and the very easy contacts 
being possible with them have created a regular exchange flow 
of requests an answers following the scheme of Pig. 2, 2nd PR, 
Chapter 3. 


Table 8 OMCr LOOK DISTRIBUTION BY CMS LAimiON TO TStXlIS TEST SITES GROUPS 














































































4. ACTIVITY PERFORMED AND RESULTS OBTAINED* 

4*1 Evaporation and soil moisture 
4*1*1 Soil covered by vegetation 

4*1 *1*2 Taylor’s expansion of the siirface energy balance 
equation (ICW) 

The surface energy balance equation vas lineari 2 ed using Taylor's 
expansion* In this vay, the original h 3 rpersurface in seven di- 
mensions could be reduced to a linear relationship betveen eva- 
poration and surface temperature, or to a surface relationship 
involving evaporation, surface temperature and albedo* 

These analytical approximations permit a rapid sensivity analysis 
of the surface energy balance equation as veil as a speedy map- 
ping of evaporation from remotely sensed surface temperatiares 
and albedo* 

Results obtained by this method from aircraft MSS data in a 
desert area vere compared vith evaporation measured on the ground. 
The standard deviation of the regression line of calculated on 
measured evaporation vas 1.0 mm/day. 

Details of the method are reported in TELLUS NEWSLETTER 15 
(Annex 3) . 


4.1, !^. Field operations 
4. 1.2.1 Flight experiment 

4. 1.2. 1.1, Germany (JRC, UG/lBW, UL/dG, UR/DG, IHI^) 

At the end of the BJMC at Pattensen (June 1''79) it vas decided 
that the various Institutes participating to the csunpaign would 
give their contribution to a report describing the type of mea- 
surements taken and the equipments used and all the technical 
details necessary to data processing (see Annex 1). 

In principle this report was thought as a guide to data retrieval 
for people participating to the campaign and it should have been 
issued by the end of 1979 without experimental data. Unfortuna- 
tely JRG did not receive the contributions in due time. 

In the meantime some contributions with experimental data were 
received and it was decided as a first step to put them in the 
report. In any case almost all the experimental data are at 
Ispra and are available to the various Institutes which partici- 
pated to the campaign. 

JRC-Ispra, which did the micrometeorological measurements on the 
test field N® 1 (bare soil) has already spoiled all the data and 
used them as input for the first run on the TELL-US model. 

Although the results appear to be reasonable they are still 
waiting to be confirmed with final data on soil thermal pro- 
perties. 

4.1 •2.2, Continuous field measurements, 

4.1.2.2.1, Italy 

Policoro campaign 1979 (JRC, UBA/iace, CSATA) • 

After the first measurement campaign in 1978 (see 2nd Progress Re- 


port, pag. 18) a 2nd measurement campaign has been carried out in 
1979 to test evapotranspiration and soil moisture models as TER- 
GRA in semiarid conditions* 

Testing the validity of such models in these particular condi- 
tions is one of TELLUS aims* 

In fact both the favourable cloud cover conditions existing in 
the Mediterranean and African areas and the more homogeneous land 
structure existing in many African regions (compared to European 
ones) make future HCMM and/or thermal satellites more suitable 
for evapotranspiration and soil moisture applications* 

This TELLUS item lines up with the Commission's initiatives in 
favour of the Associated and/or Developing Countries* It is also 
view as an initial contribution to solving overhanging deser- 
tification problems. 

Wheat was the crop choosen owing to the great importance of it. 
in April-May 1979 a measuring campaign was made on udnter wheat* 
The campaign was carried out by the JRC staff* UBA/IAGCE and 
CSATA worked under contract* 

The JRC team equipped the field of 350 :< 100 m for micrometeoro- 
Iccical and radiometer instrumentation with an automatic data 
accuisition system* 

The Acronomv Institute of the Bari University was in charce of 
the measiArement of crop-and-soil parameters* 

The CSATA, Bari had the task to staff the test area, maintain 
the equipment s , unloading and loading the magnetic recorder, and 
control that the system was working properly* The ICW, The Nether- 
lands (Dr. R*A. Peddes) acted as advisor about the preparation 
of the campaign* 10^1 is pretty interested in following data 
interpretation and in results obtained. 

The parameters measured were the same as those of the 1st measti- 
ring campaign on grass (Table 9) while some other parameters have 
been particularly stressed in this campaign (Table ID)* 


- net radiation 

- solar radiation 

- albedo 

- incident long-wave radiation (sky radiation) 

- VIS-NIR and IR radiation 

- soil temperature (at four depths) 

- heat flux into the soil (at two depths) 

- dry and wet bulb temperatures (Bowen ratio) 

~ wind, velocity and direction 

- atmospheric pressiJre 

- soil moisture (tensiometers and drilling) 

- soil characteristics (composition, pP curve, ecc) . 


Table 1C 


- rooting depth and distribution 

- crop height and crop weight vs. time 

- leaf area index 

- organic matter 


JRC’s data acquisition system proved to be highly reliable (only 
some percents of data have been lost owing to local electrical 
rains failure). 

All the data have been collected and seem to be consistent, but 
due to previous engagements (Pattensen campaign) their processing 
had to be oostponed. 
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4,2 Heat Budget 
4.2#1 Natural phenomena 

This part of the investigation concerns mesoscale heat budgets 
not particularly affected by man-made heat release* It consists 
in particular in a study of the relation between different land 
use pattern and surface temperature as well as the relation be- 
tween topography and surface temperature* 

a 

.4*2*1 *1 Analysis of thermal infrared digital data recorded over 
Belgium* (KUL/LBB) 

The data analyzed was a CCT of night pass over Belgium on May 3C, 
1978 (AA 0034*0212*0*3)* It was attempted to establish relations 
between the sxarface temperature on one hand and soil type and fo- 
rests on the other* 

Methods* The digital HCMM data was visualised in form of a com- 
puter printout* The Belgian territory represented a printout of 
2 X 1*5 m which, after photographic reduction of about 5C/i, ser- 
ved as the basic material for the interpretation* Further print- 
out maps were made by introducing class lijnits of temperature# 

By means of a zoom transfer-scope these computer printout maps 
v/ere optically superposed and compared with fS4,?ferent type of maps 
or LANDSAT images* 

The maps were road maps, soil association maps and maps of fores- 
ted areas, while the LMDSAT imagery were scenes recorded on 
17*5*76, 10*4*76, 23*3*73 and 22,5,77* The LANDSAT scenes prin- 
cipally served to identify currant land use. 

Results 

- Forests; Porestsd areas in the Ardennes could be well delinea- 
ted by their higher surface temperatures. However the difference 
between coniferous and deciduous forest types or differences be- 
tween species could not be inferred from the surface temperature 
map. (see Pig* 3-5 of Annex 2 ) 

- Soil types; A distinction could be made between v/etter soils in 

the Polder area having a high clay content and sandy soils and 
sandy loam soils. Especially dune areas showed high temperatures. 
Mere details, also in the form of maps are given in TELLUS NBV7S- 
LETTER 10 ( Annex 2 ) 


4,2#1.2 The influence oB topographic structures on niaht-^time surface 
temperatures. (UF/GI) 

The data base for this investigation vas provided by a segment of 
HCMM scene AA 0C34* 0213.0.3 of May 30, 1978* It encompasses the 
Upper Rhine Valley between Basel euid FranRfurt and the surroun- 
ding highlands* 

The segment was rectified with respect to Gauss-Krttger coordinates 
and individual sections were superposed on the corresponding sheets 
of the Official Topogi’aphic Map 1:200*000 of the Federal Republic 
of Germany (TUK 200) • Details of the procedure were described in 
TBLIiUS NEV/SLETTBR 12 (Annex 5 of 2nd Progress Report). 

Geometric Resolution and Suitable Map Scales* Under the assumption 
that the minimum diameter of map elements in thematic and topo- 
graphic maps is 0,3 mm, the geometric resolution of the VHRR 
(Very High Resolution Radiometer) of the HCMM of 600 m corresponds 
to the scale 1:2. 000. COO. For this reason, we were somewhat scep- 
tical wether an enlargement to the scale of the T&K 200 and the 
superimposition of the enlarged image with the various levels of 
information provided by this map would lead to an increase in 
knowledge. 

The material presented in this study shows that HCMM harbors an 
abundance of information which is only revealed by means of this 
kind of processing. This is especially true for the correlation 
between surface temperatures and the water network as well as of 
the relationship between surface temperat’^res and the settlement 
pattern. 

However, we should note that the resolution of the HCI^t is not 
sufficient to determine the boundaries of various geographic units 
or to establish the exact limits of topographic structures by 
means of the HCMM image. 

Nevertheless, it is essential to state that the thermal pattern 
of an HCMM image at this scale can be related to and explained 
by topographical structure. 


Night-time Surface Temperatures and Relief* A veil kno'vn rule 
of topociimatology is the thermal tripartition of hilly or moun- 
tainous areas at night; cold elevations, varm slopes and cold val- 
ley bottoms* The data allovs us to test the applicability of this 
lav to the behaviour of surface temperatures in various spatial 
dimensions • 

In general, only the difference between cold valleys and lowlands 
on. the one hand, and varm slopes and mountain tops on the other, 
could be verified* Valleys are considerably colder than slopes if 
covered by comparable vegetation* The water network which follows 
the valley systems is all but congruent with the distribution 
pattern of particularly cold surfaces* Thisi applies not only to 
the valleys in the highlands carved to a depth of hundreds of me- 
tres, but also to differences in elevation amounting to less than 
one metre on the alluvial and diluvial plains of the Upper Rhine 
■alley (cf Fig*’s 10a and 16c of Annex. 4 ), 

Nevertheless, the elevated areas proved to be as warm as the 
Slopes, in some cases even noticeably warmer (Soonwald, Pfaiser 
'./aid, Westtaunes, edge of the Black Forest near Offenburg) in 
most of the image sections presented here* Cold elevated areas 
can be found in these samples only in places where one of the 
fclloving three conditions is met; 

Grassland and cultivated fields in altitude above forested 
slopes 

Extended plateaux (cf Fig* 20a of Annex 4 ) 

Extended surfaces in areas with small differences in eleva- 
tion (cf Fig* 9b of Annex 4 ) 

At the mesoscale, if we do not consider each mountain or valley 
individually, but designate the whole Upper Rhine Valley as the 
"cold lowland", the entire, strongly dissected highland fringe 
as the "varm slope" and the highland plateau which decline slight- 
ly towards the edges of the image, as the" cold elevation", we can 
establish virtual congruence between the model and reality* 

Surface Temperatures and Forest Distribution 

It is generally agreed that forests do not cool off as rapidly 
at night as do unforested areas and that the]/, consequently, pro- 
duce higher soorface temperatures at night than adjoining grass- 
land or cultivated fields. 


Therefore, it is siArprising that a much more complex picture re-^ 
suited from this study: 

- Forests on. steep slopes or convex terrain shov very high 
temperatures (as varm as city centres)* 

- All forests on the plateau of highlands display slightly lower 
temperatures but are definitively wzamier than scattered clea- 
rings or adjoining unforested teri'ain* With the exception of 
the Upper Rhine Valley, tlie contrast between forested and un- 
forested terrain is the most esserrcial factor of the distri- 
bution of warm and cold s'orfaces everywhere, unless, of co\n:*se 
the relief dett^mines the thermal pattei\n* 

- In the Upper Rhine Val?.ey, the picUire is very complex* The 
residual forests on the holocene flood plain of the Rhine are 
depicted as being relatively warm but only in places where 
they have not lost their characteristic features due to the 
canalisation of the Rhine River and the subsequent drop of the 
ground water level* The forests on the diluvial and alltivial 
fans in front of the Pfaiser Wald and the northern Vosges 
Mountains also seem to be warmer than their surroundings* On 
the other hand, most of the Hardt forests on the Rhine’s 
quarternary accumulation plain and in part on the gravel fans 
of the rivers of the Black Forest, are as cold as, or even 
colder than, neighbouring unforested areas. This is most evi- 
dent in the southern part of the Upper Rhine Valley north of 
Basel* 

If one ignores the behaviour of the residual forest on the Rhine’s 
holocene flood plain, a general tendency for great contrast be-’ 
tween forested and \mf ores ted areas on convex landforms and for 
reduction of this contrast between the two types of surfaces in 
concave terrain prevails* As we have stated repeatedly, the highest 
surface temperatures are found on forested mountain ridges and 
steep slopes* 

There must be several reasons for the distinctive behavioi.u' of 
forested surfaces. 


- The low surface temperature of the Hardt forests can only 
be understood if one assvunes that there is a relatively homo- 
geneous layer of cold air extending beyond the actual height 
of the trees (inversion in the boimdary layer)* 

- The differing behaviour of forests on plateaux can then be 
explained by the fact that a sufficiently extended layer of 
inversion could not be created there because of stronger dyss. 
namic turbulences caused by winds at relatively high speeds* 
The sensible heat flux to the forested surfaces thus remains 
greater than on adjoining grassland and arable fields* 

- All this signifies that the much higher surface temperatures 
of forests at night cannot be attributed to, or at least not 
completely explained, by the capacity for heat storage within 
the forest but rather also by the sensible heat flux between 
the air and the forest* 

it 

- The e:<planation for the relatively high temperature of the 
residual forests on the holocene flood plain of the Rhine can 
be derived from the fact that, on the one hand, a greater heat 
capacity is available below the surface of the forest because 
of numerous residual water bodies and that, on the other 
hand, inversion in the boundary layer is interrupted or at 
least weakened, as in large cities. 

Cold Air Reservoirs and Cold Air Streams. With reference to the 
problem of cold air accumulation in concave areas and its outflow 
from valleys, the following statements can be made; 

- The fact that broad and shallow valley expanses are generally 
depicted as cold areas independent of their surface coverage, 
supports the theory that the thermal image can identify areas 
in which cold air collects and either comes to rest or flows 
onward very slowly down the valley. By contrast, in narrow 
valley areas, those sections where cold air escapes quickly 
are registered as warm areas* Furthermore, in contrast to 
cold air reservoirs, the areas in front of the mouths of val- 
leys were nowhere registered as particularly cold areas. 




T>: ; signif ies that basins accumulating cold air are depic- 
tea oy the satellite thermal image, but that the route taken 
by the cold air and the extension of cold air at the mouth 
of a valley is ignored* 

- It is generally agreed that elevated grassland is the main 
producer of cold air and that it is prerequisite for the 
night-time streams of cold air at the mouths of valleys ,as 
grassland gets colder than adjoining forests* Our examples 
shov, however, that extensively or completely forested areas 
can also create cold air accumulations* 

Basically, this was to expected* Prom the point of view of 
the heat budget, elevated forest areas with their higher sur- 
face temperatures must surely radiate more energy and thus 
the net radiation is even more negative than over grasslands* 
Therefore, forests have to withdraw a larger quantity of ener- 
gy from the air via the sensible heat flux and so contribute 
substantially to the cooling of the air* Forests do not pro- 
duce air that is as cold as that of grassland in corresponding 
locations but they produce it in greater quantities. 

Comparison of Airborne and Satellite Thermal Images* There are 
several advantages for working with satellite thermal images as 
opposed to aerial IR-registrations* It is much easier to rectify 
v" satellite image geometrically and to superimpose other levels 
of digitised information on it. 

An advantage for satellite thermal images is essentially based 
on the possibility of recording large areas almost simultaneous- 
ly, The recording of an area of the sise of the Upper Rhine Valley 
would require hours using the planes currently available for ci- 
vilian research and would result in data loaded with discrepan- 
cies caused by the time lag between the beginning and the end of 
the registration. Of course, if we are faced with the problem of 
clarifvina the boundaries of areas with varvina thermal behavi- 
cvo' with a decree of accuracy of less than 600 m, the actual sa- 
tellite data cannot be used. Satellite data supplied the same in- 
formation as aerial IR recistrations with correspondinc averacinc 
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for all studies requiring a survey of the thermal pattern within 
an area measuring 10 km x 10 km ore more (cf Pig. 15 of Annex 17) 
provided that sufficiently precise control points can be established 


for the purpose of geometric rectification in the surroundings of 
the area under discussion. 

Generally, this does not present any problems in areas with a 
strong relief. Satellite thermal data are much more comprehensi- 
ve than aircraft data for studies on a regional, rather than a lo- 
cal scale, since air-borne thermal images often obscure the basic 
correlation in thermal patterns because of a variety of irrele- 
vant topographical detail. 

Satellite thermal images certainly demonstrate the dependence of 
surface temperature on relief much more clearly than comparable 
air-borne thermal images. 

Further details of this study are given in TELLUS NEWSLETTER 17 
(Annex 4) . 


4*2.2 Anthropo<jeni*c heat release (UP/GI) 

The data base for this investigation vas provided by a segment of 
HCMM scene AA 0034# 0213*0.3 of May 30, 1978* It encompasses the 
Upper Rhine Valley between Basel (47®33'N, 7®36’B) and Franfurt 
(50®7 *N,8®36'B) and the surrounding highlands* 

The geometric corrections of the scene relative to the Gauss- 
Krttger coordinate system at a scale 1:200*000 were carried out 
as described in TELLUS NEWSLETTER 12. 

Surface Temperatures of Settlements* It was to be expected 
that heat islands from large cities would be depicted in the 
HCMM image because of their high surface temperature* However, 
it became evident that much smaller settlements can evoke a signal, 
under favourable circumstances, down to a siae of 2*000 inhabi- 
tants. 

Prerequisites for this phenomenon cire a densely built-up village 
centre and placement in a horiaontal, unforested environment* 
Otherwise, the weak signal emitted by the settlement would be 
masked by other effects on the thermal pattern* Bspecially the 
thermal effect of such small villages situated in concave terrain 
is cancelled out by the influence of the relief* 

However, not nearly all of the villages with the above mentioned 
favourable conditions are represented by a warm spot on the HCMM 
image, as the recorded phenomena are very small with reference to 
the means provided by geometric as well as radiometric resolution. 
The effect is of the same magnitude as the random noises produced 
by VHRR. Thus it is possible that the heat island of a village 
can be simulated or suppressed in each individual case* At the 
same time, the identification of a . village depends on wether its 
area is depicted in a single pixel or wether it is distributed 
on several picture elements, A signal on the thermal image can 
be expected only in the first case* 


The situation is less equivocal when we consider larger villages, 
and smaller towns with an average size of approximately lO.OCCin- 
habitants. These are clearly identified as heat islands by reason 
of their elevated surface temperature* The core area of every 
larger city with tens of thousand of inhabitants is registered 
as a heat island with temperatures a few degrees higher than 
its s-urroundings • 

The difference in temperature between the centres of cities such 
as Mannheim/Ludwigshafen (cfr.Fig.ia of Annex 4 ) and Strasbourg 
(cfr.Fig.13 of Annex 4 ) and their environments increases to 

over 4®C. 

The area of higher svu?face temperature extends beyond the urban 
built-up area to the surroundings because of the flow of warm 
air in the boundary layer. The suburban area is much warmer on the 
leeward than on the windward side of the cities* 

Industrial Complexes* On comparing different industrial areas, 
it becomes apparent that different levels of energy consi-miption 
are identified on the thermal image even though it is impossible 
to quantify the respective heat emission of an automobile plant 
and a chemical factory from the difference in the temperature 
reading. It is difficult to interpret these readings definitively 
because the thermal image by itself cannot be used to establish 
decisively wether the high temperatures of chemical plants regis- 
tered by the HCI'IM are really higher surface temperatures or 
wether they mtist be attributed to the masking effect of warmer 
exhaust gases containing aerosols. 

For further details of this study you are referred to TELLUS 
NEV7SLETTSR 1? (Annex 4 ). 


31 


e 


4*3' Satellite data calibration and atmosphere corrections . ^ 

4.3*1 Validation of HCMR Calibration 
4.3*1 .1 Vater surfaces 
4. 3*1. 1.1 Lake Geneva (JkC) 

Lake Geneva is one of the largest lakes of Surope with a surface 
of 581 bn and a maximum depth of 310 m. The availability of 
water surface temperatures and a considerable number of suitable 
HCMI-l scenes make it a rather attractive reference target. 

Ground reference temperatures. The International Commission for 
the Protection of Lake Geneva measures lake temperatures at month- 
ly intervals at 15 points in the lake. Surface temperatures from 
1C points in the “Grand Lac* were used in this investigation. 

The mean lake temperature was calculated for each month and a 
seasonal mean temperature C’orve was constructed by linear inter- 
Dolation. This c’urve, and the curve of + one standard deviation 
are shown in Fig. 2 and supplied the ground reference temperatures. 

KC>!P: scenes. CCT's of eight satellite scenes between June 3 and 
Oct. 11,1973 were employed. The lake surface was delimited on a 
video screen and the mean surface temperature and its standard 
deviation was calculated. In most cases, the thermal image could 
directly be used for this purpose, except for autumn scenes where 
contrast between water and land was not sufficient. The visible 
channel was employed in these cases. 

Slight problems were posed by the contrails of aircraft as a 
major European N-S air lane passes over the lake. Contrails as 
well as clouds were excluded from the temperature calculations. 

The scenes employed are listed below; 


Date 

Scene 

3.6.78 

0038.0148.0.3 

17.7.78 

AA 0082.0208.0.3 


.AA 0082. 1305. 0.1 ,2 

28.7.78 

kk 0093.0210.0.3 


A-A 0093.1307.0.1,2 

22.8.78 

AA 0118.1232.0.1,2 

18.9.78 

.AA 01^5,1237.0.1,2 

Q 

• 

• 

CO 

AA 0198.1221 .0.1,2 


TEMPERATURE. 


lO.OH 



5.0-^ 


QQj 

100 


JULIAN DATE (1978) 


Fig. 2 : Seasonal trend of the sxirface temperature of 
Lake Geneva. The stippeiad lines and the bacs 
indicate the standard deviation. 




Atmospheric corrections* ** These corrections v/ere calculated v/ith 
the KADTKA model using a coefficient K2 of 10* Atmospheric soun- 
dings from Payerne, 30 km of the Lake, were taken from the 

European Meteorological Bulletin and served for the calculations. 
Only the CC GMT soiandings are published in the Bulletin and these 
had to be employed to correct both day and night temperatures. 

The 12,00 GMT air temperature at the surface vas inserted for 
the day time corrections, 

Smissivity, The radiometric temperatures v/ere converted to phy- 
sical temperatures by assuming an emissivity of water equal to 
0.9815 (Becker et al, 1980), 

Results, Satellite and ground measured lake temperatures, as well 
as the atmospheric correction are presented in the table below, 
while Pig, 3 gives the seasonal trend of the two temper atiu*ss 
together ’vith their standard deviation. 

TABLE 11 


DATS 

TIME 

KCM-I* 

ATM. CORE. 

«c 

KCMI-1 

oc»-» 

GROUND 

3.6.78 

night 

13.0 

+ 1 .0 

15.7 + 1.0 

13.8 + 1.2 

17.7.78 

night 

12. C 

- 0.5 

12,8 + 0.7 

19.7 + 2.0 


day 

14.1 

- 0.1 

15.4 + 1.5 


28.7.78 

night 

11 .0 

- 4.1 

8.2 + 0.7 

20.0 +1.6 


day 

14.4 

- 2,3 

13.5 + 1.4 


22.8.78 

day 

15.7 

+ 0,3 

17.4 + 1 .2 

20.1 •+ 1.2 

18.9.78 

day 

11.8 

- 1.3 

11.8+ 0.7 

17.6 + 0.6 

iq,'ii,78 

day 

3.6 

0.0 

• 

00 

o 

• 

11 .3 + 0.6 
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* radiative temperature 

** temperature after correction for emissivity and atmospheric 
effects 





Sxcept for June 5, satellite temperatures are lower than Ground 
temperatures, while the standard deviation from the mean is ap- 
proximately the same for both, or somewhat lower for the satel- 
lite temperatures* 

The standard deviation for night data is lower than for day data* 
As effects of the scan angle are negligible at the dimensions 
of the lake, the standard deviation for both the satellite 
and ground data is a measure of spatial variability* 

In Pig* 3 , the difference between ground and satellite measured 
lake surface temperatures AT , is plotted as a function of days 
after launch* 

On the basis of this data it is certainly difficult to reach a 
definite conclusion on the behaviour of the HCM Radiometer. 

Only one day point, 38 days after launch, corresponds with the 
corrected calibration which was applied on the basis of the Vhite 
Sands data acquired somewhat earlier b'* NASA. 

The remaining five^T (day) range between 2*7 and 6*5®C and, in 
view of the rather important standard deviations, one could as 
well deduce a constant understimation of the time temperature 
by about 5®C or an oscillation of the calibration, the assumption 
being made that the slope of the calibration curve remains con- 
stant with temperature. 

The higher values of (night) seem in part to be due to the 
lowering of the stirface temperatitre of the lake at night, but 
no measurements are available at the moment to confirm this hypo- 
thesis . (see also TELLUS Newsletter 18, Annex 6) 


4*3 #1*2 Land surfaces 


4,3.1 #2#1 Pine forest (UF/GI, atmos.corr, JRC) 

A 1 X 1 Ion pine forest in the Upper Rhine Valley near Hartheim 
(47®12*N,7®36VS) could be identified on a night thermal image 
and served for this calibration attempt, 

KCt4M scene; AA 0034.0213.0,3 of May 30, 1978 

Ground measurements; The energy balance of the pine forest was 
monitored at hourly intervals by the Meteorologic Institute of 
Freiburg University, 

As the incoming and outgoing long wave radiation was measured 
above the canopy, the s^t^face temperature of this canopy could 
be calculated, 

/vH emissivity between 0,97 and 0,98 was assumed for the pine fo- 
rest. 

Results; Interpolating linearly for the exact time of the satel- 
lite overpass a surface temperature between 9,0 and 9,2®C was 
calculated which agrees with the satellite temperatuire of 8-9®C . 

Atmospheric corrections calculated with the RADTRA model and ra- 
diosonde data from Stuttgard showed the effect of the atmos- 
phere to be less than one tenth of a degree centigrade. 

Further details are given in TELLUS Newsletter 18 (Annex 6) 


4«3*2 Atmospheric corrections 

4*3«2«1 Experience with the RADTRA model (JRC) 

Serious difficulties vere encountered in the routine use of this 
model* 

It vas noted that very slight changes in the dev point depression, 
especially at pressures belov 800 mbar, had an eccessively large 
effect on the atmospheric correction* 

Tile inability of the interpolation routine to deal vith small 
or zero gradients seems to be the source of this difficulty 
vhich is illustrated by the folloving example: 

Radiosonde data vith 15 pressure levels; only four are indi 
cated: 


mb 

T®r 

DPD,®C I 

II 

730 

277 

5*0 

5*0 

700 

276 

3.5 

3*5 

670 

273 

3.5 

3.0 

620 

272 

27.0 

27.0 


The atmospheric correction for a radiometric temperature of 
288*3®K vas + 9*8®K vith data set I and +0*2 vith data set II, 
a difference of 10*K for a change in the DPD of 0*5® C in one 
layer • 

These problems have already been communicated to the HCMM 
Project Scientist* He has recognized that it vas necessary 
to produce and distribute urgently a nev version of the 
RADTRA model incorporating the changes suggested by TBLLUS 
Co-is* 


4*4 HCMM data simulation* Usage of filtering techniques for 
scalinq-up simulation* 

It is veil knovn that the gray level of image’s pixels depends 
on the radiance coming, in a particular spectral range, 


from the ground areas corresponding to the pixels* 

If each area is homogeneous no problem exists in its inter- 
pretation? if, on the contrary, a lot of different objects 
are contained, they contribute in different manners to the 
global gray level of the pixel* 

The problem becomes more important as greater dimension pi- 
xels are involved* i*e* the problem must to be faced to un- 
derstand the informative content of the data* 

Images from HCMM are typical example of this vith their pi- 
xels 500 X 500 m dimensioned* It is then necessary to study 
the so called scale effect to understand hov detailed in- 
formations aggregate to give rise to loir or very lov reso- 
lution data as they are detected for example from a high al- 
titude platform* 

The first question to vhich this work vas devoted iras the 
simulation of HCMM data starting form better resolution 
ones* 

The methodology vhich permits the correlation of different 
resolution data is based on filtering techniques in the 
hyphotesis that all changes due to altitude increase can 
be represented by a point spread function* These aspects, 
together vith a "practical" introduction to the filtering 
technique, are particularly pointed out all along the 
vork* 

Attention is here mainly focused on the scale effect in its 
general aspects also for its importance in problems such as 
registration betveen images and so on* 

Details of the method are reported in Annex 6 * 
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AL 

ATI 

CC 

COT 
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Eluropean communities 
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GDTA 
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HCMR 
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HOM 

HD 
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ID 

IR 

IPOV 

IGN 

JPB 

JRC 

LAT 
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